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Abstract 

We study the cosmology of perturbative heterotic superstring theory during the radiation-like 
era for semi-realistic backgrounds with initial A/" = 1 supersymmetry. This analysis is valid 
for times after the Hagedorn era (or alternatively inflation era) but before the electroweak 
symmetry breaking transition. We find an attraction to a radiation-like era with the ratio 
of the supersymmetry breaking scale to temperature stabilized. This provides a dynamical 
mechanism for setting the supersymmetry breaking scale and its corresponding hierarchy 
with the Planck scale. For the internal space, we find that orbifold directions never decom- 
pactify, while toroidal directions may decompactify only when they are wrapped by certain 
geometrical fiuxes which break supersymmetry. This suggests a mechanism for generating 
spatial directions during the radiation-like era. Moreover, we show that certain moduli may 
be stabilized during the radiation-like era with masses near the supersymmetry breaking 
scale. In addition, the moduli do not dominate at late times, thus avoiding the cosmological 
moduli problem. 
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1 Introduction 



One approach to modern cosmology at the fundamental level focuses on the study of time- 
dependent backgrounds that implement a period of inflation. Such a period of accelerated 
expansion provides explanations of the homogeneity, isotropy, flatness, large size and entropy 
of the Universe [Ij . It also gives an origin of the nearly scale-invariant spectrum of primordial 
cosmological fluctuations. Originally introduced in the context of field theory, a natural goal 
is to implement this scenario in string theory by taking into account both perturbative and 
non-perturbative quantum corrections and looking for solutions which are de-Sitter like [2]. 
Often in such inflationary scenarios, one encounters the cosmological moduli problem [3], 
where at intermediate times the Universe is not thermal but rather dominated by the energy 
stored in massive moduli. Their eventual decay can lead to problems such as excess entropy 
production. 

A drawback of the inflationary paradigm is that it does not explain how to resolve the 
initial Big Bang singularity. Another approach, stringy in origin, is to use the Hagedorn 
transition [H |5|, [6] that occurs at the ultra high temperature Th as an alternative approach 
to inflation [71 [S] . However, it is again often difficult to maintain analytical control when 
trying to implement dynamically such a phase transition. In order to address this problem, 
string backgrounds have been recently constructed which do not present a breakdown of the 
canonical ensemble description, due to a duality involving the Euclidean time circle [91 [TO]. 
A study of the cosmological behavior for a special example can be found in [TT] . 

In both approaches, inflation or the alternative stringy mechanism at the Hagedorn tran- 
sition, there is a problem associated to connection of these early times to the standard matter 
dominated cosmology. For example, one needs to stabilize moduli after which one often runs 
into the cosmological moduli problem [3|. In the present work, we find a radiation-like evolu- 
tion which connects the Hagedorn era (ending at a time t^) to the electro- weak era (starting 
at a time tw), where the standard model particles gain their masses. When considering 
this intermediate epoch {Ie < t < tw)-, it is possible to parameterize our ignorance of the 
earlier evolution by considering arbitrary initial boundary conditions (IBC) at Ie- From the 
outset at t^;, in a perturbative approach, one considers a string background defined by a 
two-dimensional conformal field theory and computes the quantum and thermodynamical 
properties of a space filling thermalized gas of string states p^. All statistical properties at 
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thermal equilibrium are derived from the underlying microscopic theory: First and second 
laws of thermodynamics together with state equations. The gas back-reacts on the classi- 
cally static background and in certain cases, which we characterize, leads to a well defined 
cosmological evolution. 

Following [131 ttH [IHl [IB], we implement these ideas within the framework of two- 
dimensional conformal field theories that define tree level superstring compactifications in 
four- dimensional flat space, with either A/" = 2 or A/" = 1 supersymmetry spontaneously 
broken to A/" = 0. This breaking is introduced by geometrical fluxes in the internal space 
[IZl mi [El |20l 1211 [22] . Finite temperature T is switched on by considering the Euclidean 
version of the model where the time is compactified on a circle of perimeter /3 = 2ttRq. 
Again, an appropriate geometrical flux along the compact Euclidean time is introduced to 
implement (anti-)periodic boundary conditions for (fermions) bosons. The advantage of us- 
ing geometrical fluxes, is that our analysis is exact in the string scale a'. At the one-loop 
level, the partition function (or free energy) is non-vanishing and implies a pressure and an 
energy density that induce the above mentioned cosmological evolution. A special role is 
played by the supersymmetry breaking scales; namely the temperature T and the universal 
no-scale- modulus $, which appears in all A/" = 1 effective supergravity theories and defines 
the supersymmetry breaking scale M($) [231 121] zgi'o temperature. 

The moduli that are participating in the spontaneous breaking of supersymmetry are 
running away, namely the supersymmetry breaking scale is evolving in time, while their 
ratios are stabilized [THIlS]. The analysis of the dynamics of the spectator moduli which 
are not involved in the supersymmetry breaking was initiated in [2SII2S]- There, it was found 
that for fij > T and M, the one-loop partition function (to all orders in the string length) or 
equivalently the free energy density J-" takes the following form during the intermediate era, 
i. e. between the time (the Hagedorn exit time) and the time (the electroweak phase 
transition). 



As a result, it turns out that the /i/'s are either dynamically stabilized at their self-dual 
(or enhanced symmetry) points, with masses of order the supersymmetry breaking scale, or 
frozen due to the expansion of the universe. The fact that the back-reaction of a gas of 
strings can stabilize a radius at its self dual point was initially discussed in [23 [2H] (see also 
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[T2l |29]). Here we see that once supersymmetry is broken, the stabihzation also arises from 
purely thermal/quantum effects. In particular, once the supersymmetry breaking scale is 
stabilized, as discussed below, the quantum effects will be sufficient to stabilize the radius 
with a mass around the supersymmetry breaking scale. 

For a large class of models, the time-trajectories of M{t) and T(t) are such that the ratio 
M/T is stabilized to some model dependent constant. The external space-time evolution is 
attracted to a Friedmann-Lemaitre-Robertson- Walker (FLRW) type cosmology correspond- 
ing to a "Radiation-like Dominated Solution" in four dimensions (RDS^) [30]. It satisfies 

RDS^ : M(t) oc T(t) oc ^ oc 6^"^^ oc ^ , for tE<t<tw , (1-2) 

a[t) Vt 

where a is the spatial scale factor and is the four- dimensional dilatonj^ Strictly speaking, 
the attractor (1.2) does not describe a radiation dominated universe but mimics it. This is 
due to the fact that the time-dependent state equation of the thermal gas of string modes 
does not converge to Stefan's law. It is only by taking into account the contribution of 
the coherent motion of the modulus M{t) that the total energy density and pressure in the 
Universe satisfy ptot = 3Ptot during the intermediate era. 

Along the radiation-like evolution, the supersymmetry breaking scale is not stabilized 
but rather falls in time proportional to the temperature. It is exactly this fact which allows 
us to evade the cosmological moduli problem. The masses of the moduli are related to the 
supersymmetry breaking scale and are thus also falling in time. The net effect is that the 
energy stored in the moduli dilutes faster than the thermal energy, and the evolution always 
remains radiation-like. In addition, the supersymmetry breaking scale M{t) is naturally fol- 
lowing the temperature scale T(t), thus the thermal relic density of supersymmetric particles 
is naturally low and we avoid the cosmological gravitino problem. At the end of the interme- 
diate era, t ^ tvK, it is still clearly necessary to stabilize the supersymmetry breaking scale 
and also the evolution of the moduli masses. Fortunately at least one known mechanism 
already exists, namely radiative electro-weak symmetry breaking which has been shown in 
certain cases to stabilize the supersymmetry breaking scale around the electro-weak scale 



^As a necessary and sufficient consistency requirement, we note that in this intermediate cosmological 
regime, the smallness of the space-time curvature scales — {a/ of' and the dilaton scales 0^, (/), and 
the evolving internal radii scales, (Ri/Ri)^, Rj/Ri, is guaranteed thanks to the "attractor mechanism" 
towards the RDS''. In particular, they are all decreasing at late cosmological times together with e"^ so that 
our quasi-static and perturbative approximations become better and better as time passes. 
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[3T] . This provides a natural explanation of the hierarchy between the string scale and the 
supersymmetry breaking scale and is further discussed in the conclusion. 

In Sect. [2} we first calculate the free energy density at one-loop for the simplest orbifold 
models with initial M = 2 and M = 1 supersymmetry, in order to uncover general behavior. 
We then discuss how these features apply to more realistic models, such as those of |32j . 
In Sect. |3} we study the resulting induced cosmologies. In Sect. |4| we find that toroidal- 
like internal directions with geometrical fluxes can be dynamically decompactified. This 
shows the space-time dimension can change during the intermediate cosmological era. As 
a starting point, we perform these analysis while keeping fixed the remaining spectator 
moduh /i/ at values much larger than T{tE) and M{tE)- In Sect, [sj we relax this restriction 
and analyze the dynamics of the spectator moduli in the presence of the induced thermal 
effective potential. In Sect. [6} we give a summary of our results, while in Sect. [7] we 
give our conclusions. In appendix A, we study the simpler case where supersymmetry is 
spontaneously broken by thermal effects only. There, we derive the very well known Stefan's 
law and black body physics in arbitrary dimension D from the effective supergravity theory 
of string theory models at finite temperature. 



2 Free energy density 

In this section, our aim is to derive the canonical ensemble 1-loop free energy associated to 
Euclidean string backgrounds. We first consider the simplest cases in order to understand 
the generic behavior of the free energy density. We then discuss how our results generalize 
to other backgrounds, including the semi-realistic models of [32J. For the computation to 
be sensible, we regularize the spatial volume by putting it in a "large box". By "box", it 
is understood a squared torus with radius -Rbox? while "large" means the sum on the 
Kaluza-Klein (KK) states is replaced by a continuous integral. We analyze three kinds of 
Euclidean models: In case (I), we take 

(I) : S\R^) X T3(i?box) X S\R^) X 51 X ^ , (2.1) 

£,2 
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with flux in the directions and 4. This corresponds to a background with A/" = 2 — )■ at 
finite temperature. Next, is case (II), the internal fiux wraps an orbifold direction, 

(II) : S\Ro) X T=^(i?box) X ^'(^^J""^' X . (2.2) 
Case (III) is a model with A/" = 1 — )■ 0, 

(III) : S\Ro) X T3(i?box) X , (2.3) 

with again fiux in the directions and 4. 

In these models, the FLRW scale factor a scales as -Rbox- The supersymmetry is spon- 
taneously broken at zero temperature, with mass shift scale M oc I/R4, by the presence 
of internal fiux. The computation is specialized to the intermediate cosmological era where 
T oc I/-R0 and M are much smaller than the Hagedorn temperature i.e. 

Ro and i?4 > 1 , (2.4) 

(and above the electroweak scale in realistic models). We also suppose that the remaining 
moduli associated to the internal directions (/ = 5, . . . , 9) introduce scales that are very 
large compared to T and M. Quantitatively, this means the internal radii satisfy 
1 1 

— < i?/ < i?o and — <t:Ri <t: R4, / = 5, . . . , 9 . (2.5) 

Rq /l4 

At this point, the reader interested in the induced cosmologies can switch to Sects [3] and 
|4| where the attractor solutions are found. It will be shown in Sect. |5]that the hypothesis 



(2.5) is consistent and always reached dynamically. Before we enter into our subject, let 
us comment that our formalism and ideas are presented in appendix A in much simpler 
models where supersymmetry is spontaneously broken by thermal effects only. There, the 
computation of the free energy density, the induced cosmology and the moduli stabilization 
for this simple case are presented. 

2.1 AT = 2 ^ models: Cases (I) and (II) 

In Einstein frame, the free energy density J-' is related to the Euclidean partition function 
Z as 

= " (e-*/3) (e-3<^Hox) ' ^^'^^ 
5 



where (3 = 2ttRo, Vbox = (27ri?box) and (j) is the dilaton in four dimensions. In case (I), Eq. 



(2.1), its derivation for the heterotic string can be found in fl5\ and is summarized in 



appendix B. Let us give more details in case (II), Eq. (2.2). 

To compute the Euclidean partition function, the temperature is introduced by imposing 
non-trivial boundary conditions. This is achieved by inserting a phase in the i) lattice of 
zero modes of S^{Ro). For the pure KK states, this phase is (—1)", where a is the fermionic 
number. Requiring modular invariance then fixes the form for general states. Similarly, for 
spontaneous supersymmetry breaking by an internal flux along S^{R4), one inserts in the 
associated i) lattice a phase determined by a choice of R-symmetry charge a + Q^. The 
resulting partition function is 



^ ^ H,G a,b ' no, mo 



3 I aTiaT2 1 \- 1 ^a+b+ab '^ibl ^Ib+Gl^ib-Gl ^-^Irho+nor]^ f-^arho+bno+rhono 



^ \^flm8 ^ nn']<inn'+H]nn'-H^ yia+Qi,b+L4)iH] 7(0,0) rHn ^(2,2) . 
7,5 7', 5' 



where Z^^'^^ is the standard contribution associated to the Z2-twist on T" given in (|A.4|) and 
the orbifold block Z'>'^'^9'^'''~^^'^\q] computed in appendix C is 



^(1,1) [g\ 




2. / .i_h!1.i_hi ^(a+Qi)G+{b+Li)H+UHGfimod2 fo^ {g\ ^ [o]- 

^[l-Gmi-G] 



(2.8) 

In this expression, a generic choice of charge Q4 is a linear sum of gauge group right moving 
lattice charges and/or orbifold twists, which can be represented as 

Qi = m + f]'i + f]"H , L4 = f]6 + r]'6' + f]"G , e4 = l-r7-r/', (2.9) 

where r/, f/', f]" are constants equal to or 1. L4 and 64 are determined by modular invariance 
by considering the transformations r — — 1/r and r — r + 1, respectively. 

Clearly, the contribution of the sector [q] = [q] is identical for the backgrounds (I) and (II) 
and its computation in the intermediate cosmological era where T, M <^ is summarized 
in appendix B. The result involves the untwisted "light" states only, as compared to Th- 
Quantitatively, the substantial contributions arise from the massless modes together with 
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their KK towers associated to the large radii Rq and i?4, as compared to the string length. 
Exponentially suppressed terms in Rq and/or R^ can be neglected. 

The situation is different for the contributions of the sectors [q] ^ [{]] . For the background 
(I), the dependence in R4 is similar to the one arising in the sector [[j] (see appendix B). This 
is due to the fact that the orbifold acts on directions orthogonal to S'^(i?4). However, this is 
not true anymore in case (II), where these sectors are independent of i?4, as can be seen in 



the second line of Eq. (2.8). The result is that these contributions are similar to the ones 
that would arise in a pure thermal case i. e. with no spontaneous breaking of supersymmetry 
arising from the internal space, and discussed in appendix A. The only difference between 



the [q] ^ [q] sector in (2.8) and its analog in the pure thermal case is the Kronecker symbol 
6 that may eliminate some sub-sectors. More specifically, for a choice of charge Q4 (2.9) 
with f] = 0, the (5-symbol cuts contributions which are already vanishing (they all contain a 
6[\] factor): We are back to a standard derivation of the pressure arising in the pure thermal 
case. On the contrary, for f] = 1, the 5-symbol eliminates the sectors [J] = [\Zg]- This has 
no consequence in the untwisted sector [q] = [5] which is exponentially small in the regime 
Rq ^ 1. However, in the twisted sectors [q] = [q], the subtraction of ^OQ]^ = ^0[IZq]^ from 
the Eg lattice contribution is equivalent to replacing its g-expansion in the pure thermal case 
as 1 + C(g) -> i + C(g). 

Altogether, the 1-loop free energy density J-" is given by 

(II) : T=-T' (n^ Mz) + n^l - I) C4 + K Mz)) , (2.10) 



where 



Ro M 1 1 , , 

— = — , T= -, M = 2.11 



71^ 



and C4 = — is defined in Eq. (A.7). {n'^j) is the number of untwisted (twisted) massless 
48 



boson-fermion pairs in the parent model, before spontaneous supersymmetry breaking and 
finite temperature are switched on. There is a generic contribution to coming from the 
states arising generically in the moduli space of x T^/Z2, while an additional one, A^enhan, 
arises at special points of enhanced symmetry (see [25] and Sect. 5). riy = YTs=i{~)^^'^^^ 
depends on the choice of (54-charge operator. For the specific model considered at hand, one 
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has 



= 4 [504 + iVenhan] , = 4 ■ 512 , = 4 248 + (-)^ + {-f 128 + iVe„han 



(2.12) 

We emphasize that the form of the free energy is the same for all choices of R-symmetry 
charge Q^, with only the constants and riy changing. The functions friz) and fv{z) arise 
from the relatively light Kaluza-Klein towers associated to the circles S^{Rq) and S^{Ri), 
and are given by, 



,42 



fy{z) = e^^ hi-z) . (2.13) 



k,M e2-(2A;o + 1)2 + (2^4)' 



5/2 



For comparison, the form of the free energy density for the background (I) as computed 
in [23] is (see appendix B), 



(I) : 7= -T' ((n^ + 4) hiz) + « + n'y) fv{z))., 



(2.14) 



where n^, n^, riy are given in Eq. (2.12) and riy = X]"=i(^)*^**'^'' takes the value 

■ (_)r4 (i + (_)^' ) 256. 



(2.15) 



2.2 A/" = 1 ^ models: Case (III) 



We can extend the previous considerations to backgrounds of type (III), Eq. (2.3). The 
first Z2 twists the directions 6,7,8,9, while the second Z2 acts on the directions 4,5,6,7. 
Again, = 1 supersymmetry is spontaneously broken by the discrete deformation imposed 
by the non-trivial boundary conditions along the directions and 4. Including the correct 
phases associated with the temperature and the spontaneous breaking of supersymmetry, 
the partition function is given by 

3 / dndr. 111 ..^.■,.. g[g]g[^:g]g[^:g]g[rgra] 

^-^oi?boxy ^2 ^ 2 ^ 2^^") ^5 



HuGi H2,G2 a,b 



ttR; 



2 



1 "^0 +"0 T P / _ \ amo +f)no +»Tio no 



7,5 7', (5' "o>mo 

7(a+Q4,fe+£4)rHi] 7(0,0) rHi] 7(0,0) r//2l ^(0,0) r-/fi-//2l ^ (cy ^ p\ 

^(1,1) biJ ^(1,1)[GiJ ^(2,2)b2J ^(2,2)i-Gi-G2J ;^2^ ' l^-^'^J 



where the charges Q4, L4 are given in (2.9) and the blocks z'^'^~^^^'''~^^'^\q] and are 



given in (2.8) and (A. 4). 



- Clearly, the contribution of the [q^] = [[]] sector is half the result found for the back- 
ground (I). 

- For the sectors [q^] ^ [[]], the arguments used for the background (II) apply identically. 
If f/ = 0, the 5-symbol in z''^'^^^*'^~^'"^\^] has no effect and can be forgotten, in which case 
the contributions are identical to the pure thermal case. U fj = 1, the only consequence of 
the (5-symbol is to eliminate the subsectors [|] = [\ZgI] from the answer found for the pure 
thermal case. 

- For [q^] ^ [q], the only subsectors that are not exponentially small (or vanishing) are 
[^^1^^] = [o|(i.J or [gJgj], which happen to be equal. Thus, both imply contributions to the 
free energy density which are half the result found from the twisted sector of the background 



(II) in Eq. (2.10) 



Altogether, the pressure of the = 1 — )■ model at finite temperature is 

J- = -T' ((4'" + 4") Mz) + (n^'* + 4*) (1 - I) C4 + (n^ + 4") M^)) , (2.17) 

where the integer coefficients with index T (or index V) count +1 (or (—1)*^'') for each 
massless boson- fermion pair in the ((un)twisted,(un)twisted) sector of the Z2 x Z2 parent 
model {i.e. before temperature and spontaneous supersymmetry breaking are turned on), 

= 2 [504 + Ar,„han] , n^'* = 2 ■ 512 , 

n^" = 2-512, n^* = 2-512, (2.18) 



U.U 



248 + ({-f + {-f) 128 + iVenhan] , 4" = (-)'"2 (l + {-)"') 256. 



In the intermediate era, the above form of the free energy density naturally generalizes 
to a large number of heterotic models with A/" = 1, 2 initial supersymmetry. For any four- 
dimensional orbifold construction, with the restriction that the supersymmetry breaking flux 
wraps only one internal direction, the free energy density takes the general form: 

J- = -T^ Mz) + % C4 + nv Mz)) , (2.19) 

where the integer coefficients satisfy 

n-r > , -riT < ny < ut , hx > ■ (2.20) 
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The origin of hx C4 is the contribution of the twisted sectors in which the internal direction 
associated to the supersymmetry breaking is also twisted by the orbifold action. Indeed, in 
these sectors, the "boson-fermion" mass splitting vanishes so that the contribution to the 
free energy is purely thermal. 

When more internal cycles are wrapped with supersymmetry breaking fluxes, the free 
energy depends on = M/T and additional complex structure moduli. For example, in 
the case when the supersymmetry breaking flux wraps two directions with radii R4 and R^, 
which are taken to be large compared to the string length in order to avoid Hagedorn-like 
divergences, the free energy yields a potential for the no-scale modulus M = \ / {2ti ^/R^R^) 
and the ratio R^^/ R^ p^. The general form of J-" involves a sum of integer contributions 
dressed by generalized functions /t(v)- In general, our approach can be applied to any string 
background based on a free field CFT such as free fermionic constructions, including the 
semi-realistic Pati-Salam models of 



We stress that (2.19) is valid only in the intermediate era, while at later times infra-red 
effects will become important. These infra-red effects are highly model dependent, while our 
conclusions in the intermediate era are robust. Namely, from this general form we shall find 
an attraction to a radiation-like era where the supersymmetry breaking scale is not constant 
but falls with a specific behavior determined by the expansion of the Universe, as well as a 
mechanism for stabilizing moduli. 

In order to appreciate better the above form of the free energy, one has to compare it with 
that obtained via an effective field theory approach (with initial M = 1 supersymmetry). 
Indeed, the result of the effective field theory approach {EFT) suffers from many ambiguities 
related to the UV cut-off scale A^.^: 



- J'eft = T^n* C4 + mm C4 In ^ + ^4 1 + C2 



CO ^ rj-i 



(2.21) 



T 

where n* is the effective number of light degrees of freedom below a given temperature T. 
The term (C4 In + ^4) is the re-normalized effective potential, (the renormalization 
point /i is taken here to be the temperature scale, /z = T). The term re- normalizes the 

logarithmic divergences, M'^C4log^, while the quadratic divergences M^A^^ are always 
present in any effective field theory with an initial supersymmetry A/" < 2, where C2 is 
proportional to StrAi"^. Thus, the field theory approach is unable to explain in simple terms 
the exact string result, even though it can be understood why n* — )■ n*{z) becomes a function 
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of = M/T once all Kaluza-Klein modes are taken into account (up to the scale Aco). In 
fact, the exact string computation used here and in previous works [HI [13 [2S| is necessary 
to take into account the gravitational, the observable and the hidden sector contributions to 
T (they all contribute to the integer coefficients such as ny, ny, ut, etc..) and to justify 
the absence of any M^M^^^.- term. 



3 Attraction to radiation-like dominated eras 

In the previous section, we first computed the free energy density J-" for the simple back- 



grounds (I), (II), (III) defined in Eqs (2.1)-(2.3), when the Universe evolves in the inter- 



mediate era. Supersymmetry breaking was implemented by a single geometrical fiux in the 



internal direction 4. In (2.19) we gave the form of the free energy density for general orbifold 
models where supersymmetry is broken by geometrical fiux in a single internal direction. 
The moduli associated to the directions 5, . . . , 9 are assumed to generate scales much larger 
than T and M, as will be justified in Sect. [5] In this section we proceed to analyze the 



cosmology induced by the general form of the free energy density (2.19). We will find an 
attraction to a radiation-like evolution, where the quantity = M/T can be stabilized. 
This provides a mechanism for generating the hierarchy between the string scale and the 
supersymmetry breaking scale. This result is robust and remains valid for the semi-realistic 
Pati-Salam models of [32]. 



The behavior naturally splits into two types: the first one is illustrated by case (I), (2.1), 



and corresponds to having nj- = in the general form given in (2.19). The second one 



is illustrated by cases (II) and (III), defined in (2.2) and (2.3), and corresponds to having 
ut > 0. For case (I), the study of the back-reaction of the free energy density on the initially 
static background was initiated in [30]. For cases (II) and (III), an important role will be 
played by the additional contribution from the twisted sectors, the net result of which is to 
strengthen the attraction to the radiation-like evolution. 

We are interested in isotropic and homogeneous cosmologies. More specifically, we take 
the dilaton (p as well as -R4 to depend only on time, while the space-time metric is assumed to 
be of the form ds"^ = —dt^ + a(t)^ [{dx^Y + {d^'^Y + [dx^Y]. The remaining scalars are fixed 
to constant values, while the gauge fields are taken to be pure gauge. In Einstein frame, the 
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four-dimensional low energy effective action takes the form 

where we have introduced the following notations 
$:=Ar(0-lni?4), 01 
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V3 



(20 + lni?4) 



(3.1) 



(3.2) 



J" 



M 

Y 



Rq 
i?4 



M 



ev 2 



T 



27T 27ri?4e-<P 27ri?o 
Note that in ( |3.1 ), the free energy density J-" is only a non-trivial source for the four- 
dimensional metric and $ (or M) which involves the dilaton and scalar lni?4. For our 
FLRW metric, the contribution of the free energy density to the stress-energy tensor takes 
the form T^'^ = diag(— p, P, P, P)^^ ■, where we have introduced the energy density p and the 
pressure P. They are found using the variational principle only, 

dP 



P = T*p{z) 



p = T 



dT 



P := T'^riz) 



where 



(3.4) 



and reproduce standard expressions derived from the axioms of thermodynamics. The equa- 
tions of motion and the expression of the conservation of the stress-energy tensor can be 
found in [30]. They are more conveniently written in terms of (In a)-derivatives, using 
y = H := Hy (for any field y{t)). The conservation of the stress-energy tensor yields 



din a 



T 



._,A(z)-z 



a 



with 



AAz) 



4r — Tz 
3(r +p) 



while the Friedmann equation takes the form 

r2 ^4 r{z) 



where 



}c = ^(^Mz)l-i' 



1 ^ 



(3.5) 



(3.6) 



3 — /C(z, z 

For the scalar fields, it is more appropriate to consider equations for z, defined in (3.3), 
(rather than $) and ■ 



r{z) 



3 - /C(2;, 1, 0_l) 
r{z) 



Mz) °z + A^z) P) + A-Az) U v;(.) = 



01 + ^ ^ ^ ^ ^ 



(3.7) 



. 3 — }C{z,z 



^Note that we are in the early history of the universe and so the decay rates of the moduli $ and (f>± are 
Planck suppressed and negligible compared to the gravitational friction coming from the expansion of the 
Universe. 
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where V{z) is defined by its 2;-derivative, 

V,{z) = r{z) - Ap{z) . 



(3. 



Note that the auxihary potential V{z) should not be confused with the thermal effective 
potential —J-'. 



The system of equations (3.7) is highly non-linear and allows drastically different behav- 
iors, depending on the IBC. However, a particularly interesting one can occur when V{z) 
admits a critical point, Vz{zc) = [3Ul[TniE3]. When such an extremum exists, the equations 

Zr. (/>! 



(3.7) admit the particular solution z 



est. From (3.5) and (3.6), this evolution 



describes effectively a radiation-like era in four dimensions similar to the one quoted in the 



introduction, (|1.2|). This can be seen by noting we have 

9 



3H' 



Cr 



where 



Cr = - p{zc 



a- 2 
that results from the proportionality of M{t), T(t) and l/a(t). 



(3.9) 



1 /4 \^ 

Mit) = Tit) X = — X e-^(^^) with ait) = -CJ , <j)At) = est. (3.10) 

a{t) \3 J 

Of course, a subsidiary condition for the solution (3.9) to exist is p{zc) > 0, as checked below. 
Note that geometrically, this solution satisfies 



Ro{t) oc i?box(^) oc Riit) , 
while the dilaton decreases, as follows from the definition of 



(3.11) 



in Eq. (3.2). Moreover, 



even if the evolution (3.9)-(3.10) mimics a radiation dominated Universe, it is not radiation 



dominated. This is clear from Eq. (3.8) that is telling us that the state equation of the 
thermal bath of string states is 

p = 4P. (3.12) 



Since the mass of the KK states in the direction 4 decreases as the temperature, the massive 
KK states never decouple and the thermal system is never dominated by radiation. As 
announced in the introduction, it is only by adding the contributions pkinc = -Pkine = ^'^ 
their thermal counterparts that the state equation of the total system thermal + kinetic is 
radiation-like, ptot = 3Ptot [11 [13 [25] • 

We now turn to the problem of determining when a critical point Zc exists, in order to 
classify the models. The shape of the auxiliary potential V{z) was studied in case (I) in 
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[30] . For a given value of the model-dependent ratio ny/rtT-, three different behaviors are 
allowecQ 



ny 1 

(la) : — < , no extremum, 

rij' 15 
1 ny 

[lb) : < — < , one extremum, 

15 ut 

Tly 

(Ic) : < — , no extremum. 



(3.13) 



In case (II) and (III), V{z) gets an additional linear contribution, —z x titC/^, that unifies 
the above ranges (a) and (6) and will play an important role in the next section, 

(Ila), (116) and (Ilia), (III&) : < , one extremum, , . 

(lie) and (IIIc) : < ny , no extremum. 

Fig. [l] shows the qualitative shapes of the auxiliary potential V{z) in all cases. Some remarks 









z 


{la) 





{lla) 
(Ilia) 







1 

Zcj Z 


(I^) 


— '■ — T 
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1 

Zcj Z 


{lib) 

{lllb) 


— '■ — T 





VI 


^ 


z 


(Ic) 







VI 




z 


(lie) 
(IIIc) 





Fi gure 1: The qualitative behavior ofV{z) defined in {3.^ depends on the value of the parameter ny / nx ■ 
The ranges [a), (6) and (c) correspond to —1 < ny /ut < —1/15, —1/15 < ny /tit < and < ny/riT < 1, 
respectively. For a pressure in case (I), an extremum exists only in the range (b). For a pressure in case 
(II) or (III), an extremum occurs in the ranges (a) and (6) i.e. ny < 0. (Drawing is with arbitrary vertical 
origin.) 



are in order: 

• In the range (a) of the parameter ny/riT, we see that no matter how small nr > may 
be for the backgrounds (II) and (III), its effect is drastic since it lifts the asymptotic fiat 
direction of V{z) present in case (I). 

■^The case ny = corresponds to a pure thermal case and is analyzed in appendix A. 



14 



In all instances where an extremum exists, the sign of p{zc) is positive and the radiation- 



like era (3.10) is valid. The stability of this solution is guaranteed by showing that small 
time-dependent perturbations always die off. In particular, we find that I — )■ as a — )■ oo, 
showing that z = zH <^ H. From this we conclude that the contribution of the kinetic 
energy of z to is negligible. This is sufficient to show that at late times, any energy 
stored in the oscillations of the modulus z dilutes faster than the thermal energy. This was 
done explicitly in [30] for pressures in case (I) and remains true in cases (II) and (III). Thus, 
the solution z = Zc, (f)± = est is an attractor, and we conclude that for generic IBC, the 



solution z{t), (f)±{t) converges to it. The evolution (3.10) is thus a critical solution. In a 

^ o 

neighborhood of {z,z,(f)j_) ~ (2;c,0,0), the linearized equation of z in (3.7) shows that the 
system behaves as a damped oscillator, so that V{z) can be interpreted as a potential in this 
regime, as suggested by Fig. [TJ 

• For completeness, we signal that in the range (c) of the parameter ny/riT i.e. ny > 0, 
the analysis |30] in case (I) applies identically to the cases (II) and (III). The Universe can be 
attracted to a "Moduli Dominated Era" of contraction. However, the range of validity of this 
behavior is restricted by the fact that the perturbation of the classically static backgrounds 
increases with time. 



4 Dynamical change of space-time dimension 

In this section, our aim is to determine the induced cosmology in case (la). For such a 
background, the shape of V{z) in Fig. [l] suggests that z{t), instead of being attracted and 
stabilized to some critical value Zc, may reach large negative values and behave as a modulus 
along a flat direction. In such a regime T ^ M, the statistical effects are expected to screen 
the pure quantum effects (present at zero temperature). In fact, for small = Rq/R^, the 



free energy density (2.19) for all backgrounds (I), (II) or (III) becomes 



where 04,5 are defined in Eq. (A. 7) and the neglected terms are of order 0(^e~'^'^^^°^^^'' ^). We 
are going to see that this expression is more naturally interpreted from a five- dimensional 
point of view, as suggested by the large hierarchy R^^ Rq. 
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To avoid confusions with fields normalized in four dimensions, we denote all quantities in 
five dimensions with primes. In particular, we reinterpret i?4 as the radius of a fourth external 
spatial direction. In total, the Universe is five-dimensional, homogeneous but anisotropic, 
with Einstein frame scale factor a' in the directions 1,2,3 and h in the direction 4. Together 
with the dilaton cj)' := v^0j_/2 and the temperature T', their definitions can be found in 



Eq. (D.4). Instead of b, it is convenient to work with the field 

b 

a' ' 



J - 



(4.2) 



We derive the equations of motion in appendix D. The conservation of the stress-energy 



tensor can be found in Eq. (D.IO), the Friedmann equation in (D.ll) and the equations for 



^ and in (D.12) and (D.13). As in Sect. [3| they are written in terms of (Ina')-derivatives. 
We solve the equations in two steps. First, we suppose R0/R4, small and neglect the linear 



and quartic monomials in the partition function or free energy density, Eq. (4.1). Under 



this approximation, the thermal gas is found to satisfy the state equation for radiation in 



five dimensions, as follows from Eq. (D.5), 

p' = 4P' = AT'^ TiT C5. 



(4.3) 



Correspondingly, the four equations of motion admit the particular solution where {z, ^, (f)'j_) 
are all constants, (2:0, ^o, 0j_o)' while the metric and temperature evolutions are that of a 
radiation dominated Universe in five dimensions. 



b{t) = a'{t) X e' 



Co 



1 



X e 



where 



12 



/5 



a 



-5(zo+Co) 



(4.4) 



In appendix E, we show that for arbitrary IBC, ^, (and thus z) are converging to at 



late times. This means that the solution (4.4) is an attractor whenever we may neglect the 



subdominant terms in (4.1). 



Second, we study the effect on the dynamics of the "residual forces" that arise from the 
subdominant terms or e^^ we have neglected in Eq. (4.1). The conclusions can be inferred 



from the shapes of the potentials V{z) in four dimensions, shown in Fig. [T] In case (la), 
V{z) is slightly increasing when <^ 1, so that the residual force is expected to push z 
towards more negative values. However, the conclusions look the opposite in all other cases. 



To be quantitative, we consider a perturbation around the solution (4.4), 

z = Zo + e(^) , ^ = ^0 + £{0 > 0± = 0±o + ^(<^'x) ' 



(4.5) 
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where \e(z)\, kcol? \^(<t>' )l ^^^e (I) (e^° in case (II) and (III)) are <^ 1. At order 



one, Eqs ( |mO| ), ( [m2| ) and ( |D33| ) become 

O 3 O OO 3 /o \ OO 3 O / A \ 

Hz) = - 4 ^(0 , ^ (0 + 2 (,^«) + j = ' ^ {<^Y) + 2 ^(<^x) = ° ' (^-^^ 
where the definition of the constant C diff'ers if = or nj- > 0, 
,4.0 9 ^ ^ 1^ .... rr^ - .^0.^4 



C = e^"° 2 - — + — in case (I) , C = e"°-— in case (II) and (III) . (4.7) 
C5 \nT 15/ c^ut 



The solutions to Eqs (4.6) are 

and some consequences follow: 

• In case (la), we have C < and, as expected, the residual force (compensated by the 
friction due to the expansion of the universe) induces a small negative velocity for z that 



decreases. The perturbation (4.5) is thus more and more negligible and the attraction to the 



solution (4.4) is justified. 

• In all other backgrounds, (16), (Ic) and cases (II) and (III), we have C > 0. The residual 
force implies that z increases and the approximation ^ 1 breaks down at some point. 
This is consistent with the fact that the system is actually going back to a regime T ^ M 
and converges to one of the attractors found in Sect. |3]in four dimensions. 

• In the latter situation, the force that is pulling the system out of the regime T ^ M is 
much stronger in cases (II) and (III), as compared to cases (16) and (Ic). This is clear from 



the definition of C in Eq. (4.7) that is linear in Rq/R^ in case (I) and quartic in case (II) 
and (III), or from the slopes of the potentials in Fig. [Tjfor <^ 1. 

Our conclusion is that in all cases where ny < 0, the cosmological evolutions are attracted 



to radiation(-like) solutions. Geometrically, they satisfy the proportionality relation (3.11). 
However, important remarks need to be pointed out. In case (la), z is not stabilized but 
only frozen at late times at a value which is not uniquely determined by ny/nT, i-e. the 
properties of the massless spectrum. In that sense, z enters dynamically into a phase where 
it behaves as a modulus. Its final value Zq depends on additional data, namely the IBC 
at the exit of the Hagedorn era, and in particular the initial kinetic energies. Thus, the 
backgrounds (la) provide a natural explanation for the dynamical appearance of a hierarchy 
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Ri ^ Ro which justifies that the radius R4 should not be treated as an internal space radius 
(that determines a supersymmetry breaking scale screened by thermal effects). Instead, the 
backgrounds (la) yield a dynamical decompactification of the direction 4 and stabilize to 
five the space-time dimension. In the other cases where ny < 0, we may also attempt to 
introduce this hierarchy by adjusting the value of Zc] however, it is only at the price of a 
very artificial choice of the model that the critical value Zc would satisfy e^" ^ 1. 



5 Internal space dynamics 

In the previous sections, we have concentrated our attention on the dynamics of the moduli 
which determine the scales of spontaneous supersymmetry breaking, namely M and T, under 
the assumption that the other scales characterized by internal moduli were much larger than 



M and T. For instance, the internal radii were taken to satisfy the inequalities (2.5). Under 
this hypothesis, we have shown that after the Hagedorn era, the 1-loop free energy density 
(or thermal effective potential) is independent of them, up to exponentially suppressed terms 



(see Eq. (1.1)). As a result, we were allowed to keep these moduli frozen and look for the 
cosmological evolutions of the remaining degrees of freedom. Our aim in this section is 
to argue that the above properties of the internal moduli are actually not restrictions but 



inevitable consequences of attraction phenomena. In particular, the constraint (2.5) will be 
satisfied at late enough times, for arbitrary IBC at the exit of the Hagedorn era. 



We suppose the constraint (2.5) is valid for all radii Rj except Rq (or R^), which is kept 
arbitrary. Thus, the dependance of the free energy density with respect to Rq (or R^) is in 
general no longer exponentially small. This implies that we need to examine cosmological 
evolutions with non-trivial dynamics for this radius. Depending on the initial data for the 
radius, we shall find that it either behaves as a modulus or is stabilized at its self-dual point 
with a mass scale of the order T ~ M. For large enough initial values, the internal cycle may 
experience a period of expansion, after sufficient time this expansion halts and the evolution 



settles into the RDS^ defined in (|L2|f1 



^The question of the stabilization of the internal space in the simpler case where supersymmetry is 
spontaneously broken by thermal effects only is considered in appendix A with similar results. 
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5.1 A/" = 2 ^ models: Cases (I) and (II) 



We have considered in Sect. 2.1 backgrounds where the direction 4, associated with super- 



symmetry breaking, is either a circle (case (I)) or part of a Z2-orbifold (case (II)), see (2.1) 



and (2.2) respectively. Similarly, we can take the internal direction 6 to be a circle or part 



of the orbifold. We thus have four cases to analyze: 

case (I.i) : S\Ro) x T3(i?box) x S^R^) x S\Re) x 



rp4 



case (I.ii) : S^{R^) x r3(i?box) x S^iR^) x x ^'^^^^ 
case(II.i) : S\Ro) x T^{R,,,^) x S\Re) x x ^^^^^^^ 

£,2 

/TT ••\ Ol/O \ T^S/n \ rTi2 ("^4) X S (Rq) X T 

case (II. n) : 5 (-Rq) x 2^ (-Rbox) x x — ^ — 



(5.1) 



To write the free energy density, we separate the untwisted and twisted sector's contributions, 
J-" = J-"" + J-"*. The first part, J-"", is common to the four classes of models and is computed 
in [25]. Its explicit form. 



(5.2) 



T^iz, T], = -T^ ( {Mz) + kriz, r] - |C|)) + (fyiz) + ky{z, r/ - |C|)) 
+ n^gT{z,V, ICI) 9v{z,r], \(\) 



is expressed in terms of the scalar fields 

^ _ M_Ro 
~ T ~ R^' 



Ra 



Re, 



(5.3) 



/t and fv are defined in Eq. (2.13), while the remaining functions are 



kv{z,Ti-\C\) = e^^ kT{-z,Ti-\C\ + z) 



2K5/2 [ 27r|m6|e''-l^l-Je2-(2A;o + 1)2 + (2^4)^ 



e2-(2^o + l)2 + (2fc4) 



5/4 



(5.4) 



27^5/2 27r(e2|fl - l)e^-\'^\ J e^^{2ko + 1)^ + (2^:4)^ 



e2-(2A;o + 1)2 + (2A;4)2 



5/4 



9v{z, r], Id) = e^' gri-z, r] + z, |C|) , 

(5.5) 

where are modified Bessel functions of the second kind. In J-"", the coefficients 



and riy are defined in Eq. (2.12). They count the massless boson- fermion pairs (in the 
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parent supersymmetric model) for generic values of Rq. However, the self-dual point Rq = 1 
is special, due to an SU{2) enhancement of the gauge symmetry. The additional massless 



states at this point are taken into account by the terms in the second line of Eq. (5.2), with 



4-2, 



(5.6) 



The role of the functions qt and gy is to interpolate between the generic and extended 
spectrums. On the contrary, the role of kx and ky is to interpolate between a four dimensional 
point of view (when is close to one) and a five dimensional one (when i?e or 1/ Rq is very 
large). 

We now move onto describing the twisted contributions J-"* for each background: 
- Case (I.i) : The result for these models is derived in [25] , 

= -T' (/^(^) + fc^(z,r/ - Id)) + n'y {Mz) + ky{z,r^- |C|))) , (5.7) 



where and Uy are defined in Eqs (2.12) and (2.15), respectively. 



- Case (I.ii) : Since the twisted sector of the partition function is independent of Rq, the 



result reported in Eq. (2.14) remains valid for arbitrary 



T' = -T' (</t(^)+</v(^)). 



(5.^ 



Case (II. i) : The result is as in the pure thermal case, Eq. (A. 23), up to the R-charge 



dependent factor as in Eq. (2.10), 



^* = -T^4 (l-^) (c4 + %)) 



where 



y = z + r]-\C\ 



(5.9) 



(5.10) 



and the function k defined in Eq. (A. 24) interpolates between the four dimensional regime 



where ( is close to zero and the five dimensional one for large \(\. 

- Case (II. ii) : The twisted contribution being independent of Rq, the result reported in 
Eq. ( 2.10[ ) remains valid for all C, 



C4. 



(5.11) 
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The free energy density T is the effective potential at finite temperature for the modulus 
C, as can be seen from its equation of motion [23] , 



C + + 



5C 



0. 



(5.12) 



For the backgrounds of type (I.i), this potential has been analyzed in detail in Ref. [25] and 
presents five distinct phases. Our aim is to generalize these results to the other cases in Eq. 



(5.1) and examine the similarities and discrepancies. 
• Phase 1: Let us define a neighborhood of - 



Iby 



ICI< 



1 



and 



(5.13) 



2i?o 2i?4 

In case (I.i), one finds ({t) can be stabilized at the origin C = 0, which is a local minimum 
of J-". An RDS attractor exists if the model-dependent integer coefficients in J-" satisfy the 



inequality (la) or (16) in Eq. (3.13). However, note that the relevant spectrum at the 
self-dual point includes the vector multiplets of the SU{2) enhanced symmetry point. 



+ + n*y 1 
la.i : — ^ < 

^ ' nVp + n)^ + n\, ~ 15 



(I6.i) : _l<!^ll±<±4:<0. (5.14) 



In the latter case, z is stabilized. In the former one, z freezes at an arbitrary large and nega- 
tive value zq and the cosmology is better interpreted in five dimensions]^ These conclusions 
remain identically true for the backgrounds of type (I.ii). This is due to the fact that the 



functions and ky that make the difference between Eqs (5.7) and (5.8) are exponentially 
suppressed when ( is in phase 1 (given the fact that Rq ^ 1). 

For the same reason, the models of type (II. i) and (II. ii) can be analyzed simultaneously 



since their only difference in Eqs (5.9 ) and (5.11 ) is the function k, which is also exponentially 
suppressed in phase 1. The analysis of [25] applies to these cases and one concludes that the 
above RDS attractor with both ( and z stabilized exists if ny + fiv < 0- This refers to the 



union of ranges of parameters named (Ila), (116) in Eq. (3.14) for the enhanced spectrum. 

Expanding the free energy to quadratic order, we find that the mass of ( goes as rriQ oc 

T oc M, where the exact coefficients may be computed in each case and is dependent on 

^We have generalized the analysis of case (16. i) that can be found in Sect. 3.1.1 of Ref. [25] to case (la.i). 
For any constant z = zq such that e'^^" <^ 1, an RDS is found as long as a "residual force" that scales as 
{Rq/Rh)'^ is neglected. Taking into account this correction, Eq. (3.14) in [53] gets an additional constant 

< 0, which is similar to Eq. (4.8), with the 



term (while C = 0). It follows that = "^^+"^^ 
same conclusions. 
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the critical value of z. Thus we may stabilize all radii moduli with masses of the order the 
supersymmetry breaking scalej^ One also has that C = C-^ ^ as the Universe expands 
so that energy stored in the oscillations of C dilutes faster than the thermal energy and the 
moduli do not dominate at late times. It is only at the exit of the intermediate era, where 
the electro-weak symmetry breaking takes place, that M and thus reach their final values 
of order the electro-weak scale. We would like to stress that in models where one stabilizes 
the moduli with a constant mass from the outset, one always reaches an era where the 
Universe is dominated by moduli and no longer thermal like [3]. For the models considered 
in this paper, and in particular due to the way we break supersymmetry, the moduli may 
dilute their kinetic energy during the intermediate era so that the Universe is never moduli 
dominated even after is fixed. 

• Phase 2: When C, is in the range 

and — — < C < In-Ro and lni?4 , (5.15) 



2i?o 2i?4 

the functions fcr, ^y, ^ and qt^ gv are exponentially small and the thermal effective potential 



of C is flat. As in all conclusions found in the previous phase 1 apply to phase 2, up 
to two differences. First, C,{t) freezes at an arbitrary value along the above flat direction. 
It is not stabilized but behaves at late times as a constant modulus. Second, since it sits 
away from the enhanced symmetry point, one has to replace + + — + and 

riy + fly + Uy ^ riy + Tly . 

• Phase 3: When the inequality 

C > ln/?o and InR^, (5.16) 

is satisfied, the functions kx, ky, k do contribute, while gx and gy are exponentially negligi- 
ble. Actually, for Rq ^ _Ro and -R4, the free energy density depends on two variables, z and 
y = z + 1] — \(\ only. It is thus more natural to study the dynamics of y instead of ( i.e. the 



relative motion of Rq compared to Rq. To be specific, the untwisted contribution (5.2) takes 
the simplified form 



4 
1^ 

^Current work is in progress to show this generahzes to ah moduH excluding the dilaton [34] 



^ = i e-y K /f (^) + K fv\^)) + e'' (n^T + O ^^^4) , (5.17) 
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where 

= "^^2 , A'\z) = e- ff\-z) , (5.18) 

while the twisted sectors yield, 

case (Li) : ^* = -T^ (e"^ « fP{z) + n'y fl,'\z)) + e^^ (n*. + O^C4 + ■ ■ ■ ) (5.19) 

case (I.ii) : = -T^ (n'^ Mz) + < fv{z)^ (5.20) 

case (Il.i) : J^* = -T^ < (^1 - |) (^e"^ cg + e=^^ + ■ ■ ■ ) (5.21) 

case (Il.ii) : J"* = -T^ (^1 - |) C4 ■ (5.22) 

To reach these expressions, we have neglected terms of order ^-'^'^(.Ro/Re) \ jj^ this regime, 
we find in all cases that the Universe is attracted back towards phase 2. That is, even though 
R^it) is increasing, it turns out that -Ro(^) and R4{t) end by increasing faster and always 



catch RQ{t) so that we enter into the regime of phase 2 as defined in (5.15). This attraction 
is a result of the sub-dominant terms in the free energy density]^ The interested reader can 
follow the detailed derivations and conclusions in each situation in Appendix F. 

• Phase 4: The dynamics in the range 

- lni?o and - lni?4 < ( < ^ and ^ , (5.23) 

2Rq zR^ 

is identical to the one of phase 2, as follows from T-duality C 

• Phase 5: Similarly, when ( satisfies 

- In i?o and - In i?4 < C, (5-24) 
which is T-dual to phase 3, the system is attracted back to phase 4. 

5.2 Af = 1 ^ models: Case (III) 

To extend the analysis of the internal space dynamics to models with TV = 1 — 0, we 



reconsider the backgrounds (2.3) of Sect. 2.2 In this case, we remind the reader that the two 



''For the backgrounds where the direction 6 is a circle i.e. for the cases (I.i) and (Il.i), in the hmit that one 
takes the initial value of Re to infinity or at least large enough so that one may neglect these sub-dominant 
terms, the evolution is attracted to an RDS in one more dimension. 
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Z2's act on the directions 6,7,8,9 and 4,5,6,7, respectively. The expression of the free energy 



density (2.17) is vahd when the inequahties (2.5) are vahd. Our aim is to generahze it when 



one of the internal radii Ri is arbitrary. Two inequivalent situations must be considered, 
corresponding to the different ways the generators of the orbifold act. In the first case, the 
generators treat both directions 4 and / symmetrically, in which case / = 5. The second 
case corresponds to when the generators distinguish the directions 4 and J, in which case we 
may take / = 6 without loss of generality]^ We label the cases as follows 

S\Ri) X S\Rq) X 

(5.25) 



case (Ill.i) : S^{Rq) x T^{R^o.) x 



case (Ill.ii) : S^{Rq) x T^R^^.) x 



Z2 X Z2 

S\Ri) X S\R5) X 



Zo X Zo 



In the spirit of the previous section, we separate the free energy density in four ((un)twisted, 
(un)twisted) sectors, J-" = J-""'" + J-"*'" + J-""'* + J-"*'*. Obviously, the {u,u) contribution is 



simply half the result found for the A/" = 2 — models, Eq. (5.2): 

^"'"(z,r^,C) = -T' (Mz) + kTiz,V- Id)) + (Mz) + ky{z,r^ - |C|)) 



+ n^'' griz^rj, |C|) + n^'" r/, |C|) 



where 



case (Ill.i) 



Rr 



case (Ill.ii) 



R. 



'5 • 



(5.26) 



(5.27) 



The generic coefficients n^'" and rit}^ are defined in Eq. (2.18), while n^'" = fij.^ = 2-2 



concern the additional contributions at the enhanced symmetry point C = 0. The other 
sectors are also easy to find: 

- In case (Ill.i), each generator of the orbifold group acts non-trivially on the direction 
6. Thus, J-"*'"" and J-""'* are independent of C, and their expressions computed for Rq close 



to one in Eq. (2.17) are valid for arbitrary C,. On the contrary, the product of the two Z2 
generators acts non-trivially on the directions 4,5,8,9. The contribution J-"*'* is then half the 



result found in the twisted sector of case (II. i), given in Eq. (5.9) 



-T^n'i' (1 



7] 



(c4 + k{y)) . 



(5.2^ 



The coefficient n^* is given in Eq. (2.18). 



®The / = 6, 7 and / = 8, 9 directions can be seen to be equivalent by noting that the product of the two 
generators of Z2 x Z2 acts on the directions 4,5,8,9. 
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- In case (Ill.ii), it is only the second generator of the orbifold that acts non-trivially on 
the direction 5. Consequently, J-""'* and J-"*'* reported in Eq. (2.17) are valid for arbitrary 

Finally, the remaining (t, u) contribution is half the result found in the twisted sector of 
case (Li), given in Eq. ( |5.7 ): 



■t,u 



-T' 4" [Mz) + kriz, V - ICD) + {fv{z) + kv{z, r, - |C|)) 



(5.29) 



where the coefficients ra^" and n^" are given in Eq. (2.18). 



We observe that the formal expression of J-" in case (Ill.ii) is identical to the one encoun- 
tered for the A/" = 2 — )• models we denoted (II. ii) in the previous subsection. We conclude 
that the discussions and conclusions concerning the dynamics of C, in phases 1 to 5 of the 
thermal effective potential are identical, i.e. the final attraction is a four-dimensional RDS 
in phase 1, 2 or 4. 

The dynamics in phases 1, 2 and 4 is common to cases (Ill.i) and (Ill.ii) i.e. the evolution 
is attracted to an RDS in four dimensions. For the former background, we only need to 
discuss the behavior in phase 3. When ^ Rq and -R4, we have 



-pt,u ^ -pu,t 



-T 



4 , ^-y 



n. 



U,U n{5) 



1 - 



C5 



-T^ [njl' friz) + n-^ fv{z) + n^'* {l - | j := -T^ k{z) , 



+ C>(e3^) + ■ ■ ■ ) 



(5.30) 

where the dots denote terms of order ^-'^'^iRo/^s) ^ -^g neglect. We observe that the free 
energy density in this regime is formally as in case (16. ii) in the previous section, up to two 
differences. First, the appearance of the constant term C5 makes the discussion simpler, i.e. 
valid for any Uy^ < 0. Second, the expression of k,{z) contains an additional term C4, which is 
positive. The conclusions are thus identical (see Appendix F): The evolution of the Universe 
is attracted towards phase 2, where we expect it to enter. Once there, it finally converges to 
an RDS in four dimensions^ 



6 Summary of results 

Working in the framework of perturbative string theory, we have studied aspects of the 

cosmology induced by finite temperature and spontaneous supersymmetry breaking by ge- 

^Our analytic study of the attraction from phase 3 to phase 2 does not apply to the situations where 
riy" + riy" < with Uy" > 0, but we expect this fact to remain true. 
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ometrical fluxes in one internal direction 4 of fieterotic models. We have considered simple 
illustrative backgrounds in four dimensions, with initial M = 1 supersymmetry, as well as 
models with initial M = 2 supersymmetry. The specific cases we considered are given in 



(2.1)-(2.3) and the generic form of the free energy density we analyzed is of the universal 
form 

J- = -T^ (ut friz) + nr C4 + ny fv{z)^ , (6.1) 

where = M/T is the ratio of the supersymmetry breaking scale to the temperature. J-" 
is parameterized by three integers: The choice of initial background determines ut and ut, 
while the details of the breaking of supersymmetry fixes ny- The models naturally divide 
into two types. Models similar to type (I), with supersymmetry breaking flux wrapping 
a toroidal-like direction, will always have ht = 0, while those similar to types (II) and 
(III), with supersymmetry breaking flux wrapping an orbifold-like direction, will always 
have fiT > 0. For models of type (II) and (III), hr is related to the twisted sector that is 
independent of R4. 

Whenever nv/nr > 0, we find that the Universe is attracted to a phase of contraction 
where our quasi-static approximation breaks down. The models of type (I) with ny/nx < 
— 1/15 induce a large hierarchy T/M ^ 1. The resulting evolution is well described by a 
purely thermal model in five dimensions. For the remaining cases, we find that z = In M/T, 



is stabilized at the unique zero of the force in (3.8), i.e. Vz{zc) = 0, and the evolution is 
attracted to an RDS in four dimensions. The latter is radiation-like in the sense that one 
obtains the state equation, ptot = 3Ptotj only after one includes the contribution from the 
motion of M{t). We summarize the results as 



(6.2) 





ny ^ 1 
riT ~ 15 


1 Tly 

15 Ut 


^>0 


(/) 


RDS^ 


RDS^ 


contraction 


(//),(///) 


RDS^ 


RDS^ 


contraction 



In general, depending on the low-energy particle spectrum and the initial space-time di- 
mension (at the exit of the Hagedorn or inflation era), some internal directions with super- 
symmetry breaking flux may be dynamically forced to decompactify. The analysis in [13] 
confirms this possibility on models with geometrical flux in two internal directions. It would 
be interesting to understand the full phase space for semi-realistic models that would prefer 
or at least admit four-dimensional RDS, and determine their corresponding spectrums. 
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In Sect. 5, we relaxed the frozen moduli restriction and allowed one of the radii Rj 
[I 7^ 4) of the internal space to become dynamical. In all cases, when the logarithm of this 
radius, (, is in the range defined by 



Id < lni?o and InR^, (6.3) 

it is either dynamically stabilized at its self-dual point ( = 0, with a mass of order the 
supersymmetry breaking scale, or freezes at any value of the modulus-like phase defined by 
the range 

^ and -4" < ICI < In-Ro and InR^. (6.4) 



2Rq 2i?4 

In both cases we find that the energy stored in the modulus is always diluted faster than the 
thermal energy, thus avoiding the cosmological moduli problem. When ( > InRo and lni?4. 



we find that Ro{t) and R^it) always evolve so that ({t) enters back into the range (6.3), 
where its motion halts. The attraction to the modulus phase is stronger for orbifold-like 
directions as compared to toroidal-like directions. 



7 Conclusions 



The finite temperature stringy setup naturally suggests a separation of the cosmological evo- 
lution in at least four distinct phases, according to the value of the temperature. Namely: 
(i) The very early phase, or even the "(Pre-)Big Bang phase", where the underlying string 
degrees of freedom are excited, or even strongly coupled. Perhaps string dualities can be 
applied to understand this phase and resolve the naive classical Big Bang singularity p2]. 
(a) The stringy Hagedorn phase, T ~ T^, where string oscillators and the thermal winding 
states must be properly taken into account. Both phases {i) and (ii) lead to a non-geometrical 
structure, e.g. the T-fold cosmologies studied recently in [35j. In these high temperature, 
high curvature and high string coupling regimes, the topology and dimensionality of the 
space are not well-defined concepts. Recent progress in understanding the Hagedorn phase 
has been made in [9l [101 [H] . 

{in) The third phase is the focus of this paper and has features similar to that of a radiation- 
like Freedmann cosmology [T3| HH |15l |16l |30] . Here the Universe has cooled down to tem- 
peratures far below Th and the effects of string massive states are exponentially suppressed. 
(iv) At lower temperatures, the effective field theory approach is valid. We are expecting 
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new phenomena such as the electroweak phase transition, QCD confinement and structure 
formation to take place. We also expect that in this phase, some dynamics becoming rel- 
evant at these low temperatures will stabilize the no-scale modulus $ associated to the 
supersymmetry breaking scale [23], realizing a cosmological, dynamical mechanism for the 
scale hierarchy, M^r ^ Mpianck- 

In A/" = 1 models, such a stabilization mechanism exists. At late times, one additional 
scale which enters the problem is the infrared renormalisation group invariant transmutation 
scale, Q, which is induced at the quantum level by the radiative corrections of the soft 
supersymmetry breaking terms at low energies |23l |3T]. It is possible, when T{t) < Q, for 
the radiative corrections to generate the potential for the Higgs and induce the electroweak 
phase transition, SU{2) x f/(l) — )• U{l)cm- This starts to be the case at a time tw, and for 
times afterwards, t > tw, the supersymmetry breaking scale M is stabilized at a value close 
to Q. Whether the correct Higgs potential is generated depends on the initial data; however, 
this initial value problem is avoided thanks to the attractor mechanism towards the RDS in 
earlier cosmological times. It would be interesting to realize this scenario explicitly in string 
theory, and to indeed identify models which produce the radiative symmetry breaking and 
also stabilize the supersymmetry breaking scale M. It would be especially interesting to find 
semi-realistic models which produce supersymmetry breaking scales compatible with current 
observations. 

In earlier cosmological times where M{t),T{t) ^ Q, the transmutation scale is irrelevant 
and does not modify our analysis. The results of this paper are thus valid in the intermediate 
cosmological history, Ie <t <ty^. These statements are correct if one assumes that there is 
no hidden sector gauge group G that confines at an IR renormalization group invariant scale 
Kg above Q. In string theory models with such a hidden sector, we expect the attractor 
mechanism of the intermediate era to be valid above and below A^, with threshold effects 
around Kg [36] . 

We have shown the existence of models (still consistent with the quasi-static and per- 
turbative hypothesis) which describe decompactifications of internal radii involved in the 
spontaneous breaking of supersymmetry. In these special cases, the cosmology is attracted 
to radiation-like dominated solutions in higher dimension. It would be interesting to gen- 
eralize this mechanism and try to generate the spatial directions of our Universe. We may 
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consider a scenario where three of the spatial directions form a small three-torus, T^. In 
the cases where the is wrapped by supersymmetry breaking flux, it is possible for it to 
decompactify and generate the three-dimensional space. One could try to realize this sce- 
nario explicitly within string theory and investigate the connections between the resulting 
low energy particle spectrum and the requirements to generate three-dimensional space. One 
candidate for studying this mechanism, when one direction has already decompactified, is 
the recent MSDS models [TO] . 

In scenarios where the Universe starts out very small, close to the Planck or string scale, 
and taking the results of this paper into account, we find during the radiation- like era that 
while the internal radii (not involved in the supersymmetry breaking) may initially expand, 
they are always attracted to their flat potential phase where their evolutions halt. This 
provides a natural mechanism for keeping moduli at or near the string scale and eventually 
stabilize them at enhanced symmetry points. 

We have chosen our setup so that the underlying two-dimensional conformal field the- 
ories are exactly known in string length. This restricted the supersymmetry breaking to 
occur via geometrical fluxes. It would be interesting to use string-string dualities and re- 
interpret our results as non-perturbative effects in a dual theory. This would allow us to 
understand the role of temperature in more general cases where supersymmetry is broken 
by non-perturbative effects. For example, in the type II version of our setup there is per- 
turbatively no enhancement of symmetry at the self-dual point. However, a dual type II 
description of the heterotic gauge group enhancement can be considered in terms of singu- 
larities in the internal space. For instance, a type IIA D2-brane wrapped on a vanishing 
CP^ cycle of radius dual to Rq can give rise to an SU{2) gauge theory and admits a mir- 
ror description in type IIB [37] . The equivalence between the brane- world and geometrical 
singularity pictures can be analyzed along the lines of [HH] . 
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Appendix A: Black body and Stefan's law revisited 

In this appendix, our aim is to present our approach in the simplest context. Namely, we 
consider at the classical level supersymmetric string models in /^-dimensional flat space- 
time. At finite temperature, the supersymmetry is spontaneously broken and a cosmological 
evolution is induced at the quantum level. Restricting our analysis to the dynamics that 
follows the Hagedorn era, we have T <^ T^. For simplicity, we first suppose that all other 
scales in the model are much higher than T. It follows that the only states that can be 
thermalized are massless and we recover Stefan's law and an attraction of the Universe 
towards a radiation dominated era. Then, we show the above hypothesis on the scales of the 
model is actually a consequence of the dynamics of the internal space moduli. 

From a statistical physics point of view, the system is the Universe filled with a thermal 
gas of states. To be concrete, the space is treated as a large box with periodic boundary 
conditions along its D — 1 dimensions, which is nothing but a torus T^~^. To regularize IR 
divergences, the sums over the KK states along these directions are replaced by continuous 
integrals. The supersymmetric spectrum of states which is thermalized is the whole set of 
string modes of a given model. This guaranties that the loop corrections are also free of UV 
divergencies. From a dynamical point of view, the gas exerts a force on the "walls of the 
spatial box" and, if the induced perturbation of its radius is small, a quasi-static evolution 
takes place. 

To compute the canonical ensemble free energy, we consider Euclidean backgrounds of 
the form 

S\Ro) X T^-i(i?box) X Mw-D , (A.l) 

where -Rbox is the radius of the circles of T^~^ and A^io-_d is the internal space that preserves 
at least one supersymmetry in D dimensions. Rq is the radius of the Euclidean time circle 
along which bosons and fermions have different boundary conditions, so that all supersym- 
metries are spontaneously broken by thermal effects. To be specific, we consider heterotic 
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(A.2) 



models with Aiio-o = ^ , whose 1-loop partition functions take the form 

7_ D pD-l f dTidT2 1 ^a+b+ab ^ibT ^(lo, ji,io_D)r(o,i6) 

g--^l™0+n.OTpj-_'^amo+bno+mono 
no, mo 

In this expression, the F's stand for Narain lattices. Orbifold models with Aiio-o = T^~^ x 
T^/Z2 have instead, 

fal 2 /) ra+Hl n la—Hl 



'2 H,G a,b 

r(6-D,6-D) r,(Ofl)rm 1 

^4^20 ^(4,4) IgJ 2 ^-^t^J 2 

7,(5 7', 5' 

ttR^ 

g-^|mo+no-rp^_'jamo+6no+mono 

no, mo 



(A.3) 



where we define 



(n,n) J 



r]r] 



V2 (A.4) 
for [^] ^ [°]. 



Due to the phase (— )^"°, the GSO projection in the odd winding sector is reversed and 
tachyons occur when Rq reaches the Hagedorn radius, close to 1. As said before, we are 
interested in the low temperature regime (compared to Th) where 

i?o > 1 , (A.5) 

so that Z is well defined. Moreover, we suppose for the moment that the contributions to 
the masses from the internal lattice are "heavy" . For instance, the internal radii Ri satisfy 

^ < i?, < i?o • (A.6) 

Following the steps detailed in Ref. p3] (or appendix B for more sophisticated models where 
supersymmetry is spontaneously broken even at zero temperature), the partition function 



(A.2) becomes 

Z^«„.c„ where ..^^Epj^. (A-^) 

ko 
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where ht is the number of massless boson-fermion pairs in the supersymmetric parent model, 
when the temperature is not switched on. These states give rise to KK towers of modes along 
the Euclidean time circle (the discrete sum on is obtained by Poisson resummation of the 
KK momenta along S^{Ro)). In this expression, we have neglected exponentially small terms 
of order Q(^q-'^'^Ro J^iass^^ where Mass stands for a mass contribution arising from oscillators 
and/or internal lattice zero modes. The integer number depends on the specific values 



taken by the internal moduli. In the toroidal case (A. 2), on has tit = 8 (504 + A^enhan), where 
Aenhan stauds for the additional contribution that arises at an enhanced symmetry point. In 



the orbifold models ( A.3 ), n-r refers to untwisted and twisted modes, = 4 (504 + Aenhan + 
512). 

In the regime we study, the notion of space-time is well defined in field theory (the radius 
of curvature is large in all D dimensions) and the dynamics at low energy can be described 
by an effective action S in Lorentzian time. The first non-trivial contribution to the vacuum 
energy arises at genus one. Supposing the string coupling e"^ in D dimensions is small enough 
for perturbation theory to be valid, we can write S at 1-loop order. 

Rat _ , „ , N o \ Z 



S 



St 



-2(l> 



+ 2{d<pY + ■■■ + 



(A. 



/5Hox. 

for 

In this expression, we 



for the string frame metric Qstfiv and the dilaton. The dots stand for the other massless 
degrees of freedom, while the massive states are integrated out, 
denote /3 = 2tiRq and Vbox = (2vri?box)^~^- For -D > 3, the action can be rewritten in the 
Einstein frame as, 



S 



d^Xy/^ 



where we have defined 

J" = -T^ riT CD , 



_ 20 5 

2-nRne 0-2 
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(A.9) 



(A.IO) 



Supposing the back-reaction of the 1-loop source T on the classical space-time induces a 
quasi-static evolution, we look for homogeneous and isotropic extrema to 5, with metric and 
dilaton ansatz, 

ds'' = -N(t?dt^ + a(t?\ idx'? + ■■■ + (dx^'-^Y I , 0± (t) 



-Nitfdt^ + a{tf [{dx^ + ■■■ + {dx^'-y] , 

20 1 20 

where N{t) = 2'kRo e 0-2 = __ ^ a{t) = 27r/?box e ^-2 



(A.ll) 



^''The 1-loop corrections to the kinetic terms can be absorbed by wave function renormalization. They 
would translate into corrections to the vacuum energy at second order only. 



32 



and trivial background for the other massless fields. Note that the laps function is by con- 
struction the inverse temperature, since it is obtained by analytic continuation from the Eu- 
clidean background. The components of the stress-energy tensor Ttot^i" = diag(— ptot, -Ptot, -Ptot; Ptot)^J.'' 
satisfy 



^tot = ?0l-^, Ptot = ^^l+J'-T% =^ p^^^ = _p^^^ + T^. (A.12) 



Separating the classical and 1-loop contributions, ptot = P± + P, Ptot = P± + P, the system 
(p, P) satisfies Stefan's law, the state equation for radiation in D dimensions, 

p= {D -1)P = {D -l)T^nTCD. (A.13) 

Note that for more general expressions of J-", we recover from the variational principle 
the usual quantum statistical results 

where Fgt is the free energy associated to the thermal partition function Z = Ti e~^^ , 

Redefining a more conventional time variable corresponding to a constant laps function 
equal to 1, the Friedmann and (f)± equations of motion are, 

l{D-l){D-2)H' = l<Pl + p, {^^^y (A.16) 
4>± + iD-l)H^^ = =^ 0^ = ^, (A.17) 

where c± is a constant. The conservation of the stress-energy tensor can be used to relate 
the temperature to the scale factor, 

p + {D -l)H{p + P) = =^ aT = aoTo, (A.18) 



where cto^o is a constant. The Friedmann equation (A.16) can be used to show that a{t) — )■ 



+00 as t — )■ +00. This implies that for late times, the classical kinetic energy density 
P± IN is negligible, compared to the thermal one p oc Therefore, for arbitrary 
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IBC at the exit time t^; of the Hagedorn era, the cosmological evolution is attracted to a 
radiation dominated era, since 

1{D-1){D-2)H^ = ^ where = P - l)(aoTo)^nTCB , (A.19) 

and the fields satisfy 

Stabilization of the internal space 



At this stage, we have supposed the space A^io-z? in (A.l) is static and all scales determined 



by the internal moduli are far above the temperature. We want here to examine if these 



conditions are reached dynamically. For concreteness, we consider the orbifold model (A. 3) 
with D = 4 and internal space Aie = x T^/7j2. Our aim is to study the evolution of 
the moduli of this space where, for simplicity, we allow only one radius of either or the 
orbifold factor to be arbitrary and dynamical. In both cases, we denote this radius as Rq, 
while keeping the other moduli frozen and such that the associated scales they define are 
much larger than T. Precisely, the internal space is taken to be either 

Me = S\R,)xS'x'^ or J.2 ^S\R^xT^^ 



and for instance the i^/'s in all internal directions but 6 satisfy the inequality (A.6). For 
convenience, we introduce the notations 



e 



^ = i?o , e"^ = i?6 • (A.22) 



In the case i?6 is a radius of T^, we use the general analysis of appendix A in Ref. [25] 
to write the free energy density for arbitrary as, 

T = -T' (k + O (C4 + k{T - Id)) + nV, g{T, |C|)) . (A.23) 

In this expression, = 4 [504 + Ncnhan] and nip = 4-512 are the numbers of massless boson- 
fermion pairs in the untwisted and twisted sectors (at zero temperature), when Rq takes a 
generic valuej^ However, an SU{2) enhancement of the gauge symmetry arises at Rq — 1. 

^^The integer -/Venhan refers to the additional massless contributions that may occur at specific points of 
the space of internal moduli other than R^. 
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The additional massless states that arise at this point are taken into account by the term 



proportional to ra^ = 4 ■ 2. In Eq. (A. 23), the functions we have introduced are 



me fco ^ ° ^j^ 24) 

where K2 are modified Bessel functions of the second kind. They depend on |C| only, as a 
consequence of the T-duality -Re — 1/-R6 of the model. 

When Rq is one of the radii of T^/Z2, the untwisted contribution is as in Eq. (A.23). 



Since the twisted sector is independent of -Re, the contribution to the free energy density 



can be computed for -Re satisfying Eq. (A. 6). The net result is 



J- = -T^ (n^ (c4 + k{T - Id)) + h^g{T, Id) + 4C4) . (A.25) 

To discuss the dynamics of R^it) in both cases, we may consider a general free energy 
density of the form 

J- = -T' l^nr (c4 + k{T - Id)) + nrgir, |d) + ut C4) , (A.26) 

where the main difference is that = if -Re is a radius of and > if it is a radius 
of the orbifold factor. Our interest in J-" is justified by the fact that it plays the role of an 
effective potential at finite temperature for -Re, as follows from the equation of motion for ( 
(see appendix B of [25]). 

C + 3HC+ — = 0. (A.27) 

As shown on Fig. [2| this potential presents five phases. They can be analyzed along the 
lines of 



Phase 1: In the neighborhood of C = defined by |d < ~f^j we have = in Eq. 

2-Ro 



(A.26), up to exponentially suppressed terms, thanks to Rq ^ 1. We also have g{T,0) = C4 



and J-" admits a local minimum for (, since 

jr = (^^riT + hT + hr) C4 - ut e'" C' + C(C^)) • (A.28) 



It follows that Eq. (A.27) takes the form of a damped harmonic oscillator, with time depen- 



dent coefficients, and qualitatively one may expect that ( is attracted to the solution C = 
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\nRn InRf 




Fi gure 2: Thermal effective potential for the scalar ( — IhRq, at t — IhRq and temperature T constant. It 
is nothing but the free energy density J- . 

for an expanding Universe. To make the argument more precise, we may consider a general 
perturbation around the point C = along the lines of [SI]. Under such a perturbation, 
we find the evolution is attracted back to the point C = 0? while the cosmological evolution 
is "radiation-like". What is meant by this is that the scale factor evolution is that of a 
radiation dominated Universe but the total energy density contains contributions from both 
thermal radiation and from the moduli. The fluctuation in sources the dilaton, resulting in 
a logarithmic behavior towards weak coupling and so 0_l does not asymptote to a constant 
finite value. Note that in higher dimensions, the Universe is radiation dominated and that 
(/)_L asymptotes to a constant finite value. 

• Phase 2: When — — < C < r, both functions h and g in Eq. (A. 26) are exponentially 

2-Ro 1 I 

suppressed and we have 



(A.29) 



Thus, any constant C,{t) in the above range solves Eq. (A. 27) and one finds that the Universe 
can again be attracted to a radiation dominated era. However, is not stabilized as before. 
It behaves as a modulus frozen by the friction term arising from the expansion of the Universe 



in Eq. (A. 27). Since RQ{t) is freezing while -Ro(''^) increases, the inequalities (A. 6) are better 
and better satisfied, even for / = 6. 



Phase 3: When C, > the function g in Eq. (A. 26) is exponentially suppressed and 



it is more convenient to study the dynamics of the field 



y = r-\C\ 



(A.30) 
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instead of (. Using the identity 

5/2 



C4 + k{r - Id) = h{y) := Yl 7 ' (^-31) 

k„rn, [e^y{2~ko + iy + ml 



which is obtained by Poisson resummation on mg in the definition of k, we have for y < 

jr=-T^(nTh{y)+hTC^) ■.= ~T^p{y). (A.32) 
The coupled equations of motion for the temperature, scale factor, dilaton and y are identical 



to the one considered in the problem solved in Sect, p^and given in Eqs (3.5)-(3.7) (together 



with the definitions (3.4) and (3.8)) under the replacement z ^ y. In particular, the field y 



has a potential whose derivative 



5r(|) ^ m2g% 



Vy = -TIT y ^ - UT C4 (A.33) 



is negative for any y < 0. We find that there is always a force pushing y towards larger 
values, meaning that -Ro(^) catches Re{t) when the Universe expands. In other words, the 
system is attracted to phase 2, where y > and ({t) starts freezing. We justify this in two 
steps. First, note that in the limit <^ 1, one has Vy = — e^^ 32c4/15 — nTC4. This 
shows that the force towards the phase y > is essentially constant when ht > and much 
milder when ut = 0- More generally, for arbitrary IBC at t^; such that y < (but not 
necessary ^ 1), we have completed our analytic study by a numerical simulation that 
confirms that y is always driven to phase 2. However, when ut = and one chooses the 
initial value of Rq to be infinite or sufficiently large, we can neglect the effect of this force. 
Then, the dynamics is better understood in higher dimension, where ^ 1 behaves as 
a modulus frozen by the friction term arising from the expansion of a radiation dominated 
five-dimensional Universe. -Ro(^) and R^it) are then running away, proportionally to -Rbox(^)- 
The proof of these statements is identical to the one given for the problem treated in Sect. 
|4| Eqs (4.1)-(4.4), under the replacement i?4 — j- Rq, z y, ny ^ ht- 



• Phase 4: The dynamics in the regime —t<C,< — is common to phase 2, as 

2-Ro 

follows from T-duality C — ^ — C- 

• Phase 5: The dynamics for < — r is T-dual to that of phase 3 and the system is 
attracted back to phase 4. 
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Appendix B 



Here, we give a summary of the derivation of the free energy density in models whose 



Euchdean heterotic background is given in (2.1) and considered in [Ml |25]. The partition 
function is 



•J F /To u ri „ u 



77* 

2 H,G a,b ' 

7,<5 7',5' ' ' 



''"^|^flo+nor|2^^^^amo+6no+^flono g- ^ |m4+n4rp (a+Q4)m4+(fe+-L4)n4+e4m4n4 



where Q4,L4,e4 are defined in Eq. (|2.9l) and Z^°n°'' is given in Eq. (|A.4l). Redefining 



rii = 2li + hi, rrii = 2ki + gi, {hi, gi = 0,1 for i = 0, 4) and a = d + ho + h^, b = b + go + g^, 
the use of Jacobi identity yields 



n\l+hi)+h4~i2n\l-ho-h4-H~in\l-ho-h4+H-i 

'^[l+go+g4J '^[l-go-g4-G J'^[l-go-g4+G J j^g 2) 

^0,90,^4,94 ^ 



E e ^'^ 

io,fco,U,fc4 



^?l(2fc.+gi) + (2/,+h,)T|2 . . . /l.^^.+ZlOgg+g^/l^g^+Q^^ 



In the intermediate cosmological era, we have Ri ^ 1 {i = 0,4) and there is no Hagedorn- 
like singularity. Thus, the contributions with non-trivial winding numbers in the directions 
and 4 are exponentially suppressed i.e. we can keep the sectors with hi = li = only. 
Among them, the sub-sectors with ^0 + ^4 = are supersymmetric and do not contribute (as 
seen from the presence of 9[\] functions). Due to factors of order e"'^'^?/^^ in the integrand, 
the substantial contributions in the r2-integral arise for T2 — ?■ +00, up to exponentially 
suppressed terms in This means that the integration over the fundamental domain can 
be replaced by the sum over the entire upper half strip. Altogether, one obtains 

J-1/2 Jo ZT2 H,G 7,5 7',5' ' ' 

^[o?^[^-gW^g] ^ (_)Q4g4 ^ g-:^[i?g{2fco+go)2+R|{2fc4+g4)2] j-g 



gO+g4=l fco,A:4 
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The last sum over ko, arises (by Poisson resummation) from the KK towers of states 
associated to the directions and 4. Expanding the rest of the integrand in series of q^q^, the 
Ti-integral implements the level matching condition. The change of variable T2 = x7r[i?Q(2A;o+ 
go)'^ + Rl{2k4 + g^y] shows that the monomials 6-2^^-2(^+5) ^j^]-^ j^qj^ vanishing A + B are 
exponentially small, compared to the massless contributions with A + B = 0. Note that in 
the untwisted sector H = 0, there are two kinds of such states. Some are present at any point 
in the moduli space of 5*^ x T^/Z2, while additional ones arise if we sit at some enhanced 
symmetry point, where states with non-trivial winding modes in the directions / = 5, ■ ■ ■ ,9 



are exceptionally massless. The final result for the free energy density Eq. (2.6) measured 



in Einstein frame is given in Eq. (2.14). 



Appendix C 



Let us derive the orbifold block ^(j^+«4.fe+i4) ^^^^ ^^.S). The form of z^^^^^^^'^^^'^ [O] 

can be understood from its analogue in a pure KK field theory. In this context, any field is 
multiplied by a phase (— )('^+Q4)m.4 -^y^gn — )■ x^ + 27ri?4m4. In string theory, each KK mode 
is the lowest of a tower of winding states and the generalized phase is found by imposing 
modular invariance, 



7(a+Q4,6+L4) rQ-i 

"(1,1) ^oJ 



1 Ra 



VV V^2 



E 



- |rn4+n4T|^ 



(-) 



(a+Q4)rh4+{b+L4)n4+e4m,4n4 



n4m4 



rjrj 



Tr 



Pl,r 



— 



a+Q4+e4n4 



Ra 



n4,m4 



12 

PLtotqZ'PRtot 



(C.l) 

T u^Ra (C.2) 
(C.3) 



The second line (C.2) is obtained by Poisson resummation and involves shifted momenta. 



while in (C.3), pLtot and pRtot include the contributions of the bosonic oscillators in the 



direction 4. 

The contribution of the sector [q] = [5] is found by inserting the Z2-generator g into 
the trace. Since g acts on the zero modes Pl,r (and the oscillators) as a minus sign, a 
state contributing to the trace must have zero winding and shifted momentum i.e. n4 = 0, 
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0. We thus have 
5. 



^(a+Q4,b+i4) rO 
^(1,1) ll 



a+Q4,0 mod 2 



5, 



a+Q4,0 mod 2 



oscillators 

2 



(C.4) 
(C.5) 



where A^4, A^4 are the number of left and right bosonic oscillators in the direction 4. 

Finally, the blocks associated to the twisted sectors [^] = [^^ are found by modular 



invariance, as summarized in Eq. (2.8). Note that the full string partition function can 
be found using a different reasoning. One can consider an S^{R^ (with trivial boundary 
conditions) on which one acts with two Z2-actions, the former being free and the latter not 
freely acting. Although equal, the final result appears in a different form [39] . 



Appendix D 



Our starting point is the same string theory backgrounds (I), (II) and (III) and partition 



functions that gave rise to the free energy in Eq. (2.19). However, we reconsider the effective 



action (3.1 ) in terms of fields and a free energy density defined from a five dimensional point 



of view (we denote them with "primes" to avoid any confusion with their four dimensional 
counterparts). Our goal is to find the equations of motion from this perspective. Specifically, 
we have 



S 



R' 



(D.l) 



2 2 

where is the "normalized" dilaton (p' and J-"' is the free energy divided by the volume in 
five dimensions, measured in Einstein frame, 
^' 2 _, Z 



20 



e-^27ri?o)(e-^V;'„J 



box 



(27ri?4) Vt 



box ' 



(D.2) 



J-'' takes a form analogous to Eq. (3.3), 

J-'' = —T'^p'{z) where p'{z) = p{z) 



Rq 
Ra 



r 



1 



_ 2^ 
277 itn e 3 



(D.3) 



The homogeneous but anisotropic metric ansatz in Einstein frame is now given by 
ds'^ = -N'{tfdt^ + a'{tf [{dx^f + {dx^f + (c/x3)2] + h{tf{dx^f , 01(t) , 



where N'{t) = 2tiRo e"^ = , a'{t) = 2txR^o. e 



20' 



h{t) = 2nR4 e 



2<t>' 



(D.4) 
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Varying S with respect to this metric, the 1-loop contribution to the stress-energy tensor 
defines the thermal energy density and pressures, T'J' = diag(p', P', P', P', P'4)^" , where 



BP' 

p' = r'p'{z), p^ = p' + b^ = r'{p'-p'^ 



db 



BP' 

-P' + T'^ = T"'r'iz), 



with 



(D.5) 



There are three independent Einstein's equations (for /i = z/ = 0,1,4) and one equation 
for the scalar (f)'j_. For convenience, we consider linear combinations of them and write the 
conservation of the stress-energy tensor and Friedmann's equation (we denote H' = a' /a, 
K = b/b), 



p' + {3H' + K){p' + P') + Kb 



dP' 

~db 



0, 



(D.6) 
(D.7) 



Defining = — , we have K = H' + C, and we choose the last two equations to be 
a' 

dP' 

^ + {3H' + K)^ = b — , 



db 

h'^ + {3H' + = . 



Proceeding as in 



dy 



we introduce derivatives y = — and use Eqs (D.5) to recast 

a In a' 1 r 



Eqs (D.6)-(D.9) in the form 



(D.8) 
(D.9) 



[r'^{z) - 5r'{z)] z + [p'{z) - Ar'{z) - pl^^)] ^ = 0, 

r/2 _ rpl^ ^ 



6 + 3^- 



o ' 
1 Ai2 



r'{z) 



6 + 3^-i02 



— i+p\z)i + p'^{z) = Q, 



r'{z) 



6 + 3^-^02 



— 0l+p'(^)0l = O. 



(D.IO) 
(D.ll) 

(D.12) 

(D.13) 



Appendix E 



We consider here the equations of motion (D.10)-(D.13), where the terms linear and quartic 



in = Rq/ Ri are small compared to 1 and neglected in the partition function or free energy 
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density (4.1 ). In this regime, our aim is to show that for arbitrary IBC, the two scalar fields 



^{t) and (f)'_j_{t) converge to constants for late times. 



For ^ ^ 0, we define a new function x such that = x^- Using (D.12) and the state 



equation (4.3), one finds immediately that x is a constant and we need to solve the single 
equation 

OO O / O -y^ o \ 

= 0. 



4e + e(6 + 3e-^e' 



(E.l) 



When X 7^ 0, the solution to (E.l) is 



|^|r+r_ _|_ _ r_{r_-r_,.) _|_ _ r,,. (r+ -r_ ) 

where r± = ^"^^^t^^^^ , r_ < ,^ < r+. 



(E.2) 



and Oq is an integration constant. For P' of the form (4.3), the r.h.s. of Eq. (D.8) vanishes 
and yields ^ = where 7^ is a constant, so that 



e 



(E.3) 



/ / /8 

As follows from Eq. (E.2) and shown in Fig. 3, this quantity versus (a /og) has two 




(aVao) 



Figure 3: ^ = T^I^ '^^ function of {a' /oq)'' . For > 0, the expanding solutions correspond to the 
positive branch, while for < they follow the negative one. In both cases, the scale factor tends to infinity 

o 

for late times and ^ goes to zero. (The indicated tangents can be either horizontal or vertical.) 



branches. The expanding solutions have a > and are described by the positive (negative) 
branch when > (c^ < 0). In both cases, we observe that when t increases, the scale 

o o 

factor is monotonic and diverges to +00, while (^, — )■ (0,0). 
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- When X = 0? the solution of Eq. (E.l) is 



a' J 



where s = sign ^ . 



(E.4) 



Using Eq. (E.3), one obtains 



a 



3/2 



(E.5) 



The expanding solutions have d' > and thus s = signc^. Dra wing d 'e^ versus (a'/a'o)'^''^, one 

o o 

concludes that a'{t) always diverges when t increases, and Eq. (E.4) implies — )■ (0, 0). 



Finally, for ^ = 0, ^ 0, we only need to consider Eq. (D.13) that reduces to 

1 



oo o 



40'^ + 0l(6--0'i 



0, 



and whose solution is 



eVl2 



(E.6) 



(E.7) 



where e = ±1. Solving Eq. (D.9) to find (f)'j_ 
yields 



(where c'j_ is a constant), Eq. (E.7) 



3 ' 



(E.8) 



where is the constant value of C,- The expanding solutions have 6c'j_ > and satisfy 

o o 

a' — )• +00 when t increases. In this limit, Eq. (E.7) implies — )■ (0,0). 



Appendix F 



In this Appendix, we complete the discussion in Sect. |5.1[ We consider phase 3 of the thermal 
effective potential of a radius Rq that is not participating in the breaking of supersymmetry. 
We argue that with initial data satisfying Rq > Rq and R^, the resulting dynamics implies the 



43 



evolution is attracted to phase 2 i.e. where Rq and R4 > Rq. The cases under consideration 



are given in (5.1). In the regime Rq ^ Rq and R4, the untwisted sector of the free energy 



density is given in Eq. (5.17), while the twisted sectors can be found in Eqs (5.19)-(5.22). 
• The case (I.i) is studied in [22], when the spectrum satisfies 

(I6.i) in 5 dimensions : < — < 0, (F.l) 

^ ' 31 n^ + n*. ' ^ ' 



which is the analogue of the second condition in Eq. (3.13) in five rather than four dimen- 
sions. Actually, it is found in Sect. 3.3 of [25] that when one neglects the subdominant 



terms e^^ in Eqs. (5.17) and (5.19), an RDS attractor exists in five dimensions for arbitrary 
y = Ho such that e^^° ^ 1. The direction six is part of the space-time and R^it) runs away 
proportionally to Rhox{t)- However, defining y := i/q + e(^y), one finds by taking into account 
the correction terms 0{e^y°) that e(^y) > (and of order e"^^"). It was concluded that -Ro(^) 
(and similarly R4^{t)) are "catching" RQ{t). We expect this phenomenon is valid until the 
Universe is attracted back into phase 2 (it would be interesting to confirm this fact by a 
numerical analysis as in the pure thermal case in appendix A). Once in phase 2, C(t) freezes. 



the final space-time dimension is four, and the inequalities (2.5) end by being better and 
better satisfied for all radii, including Rq. 



This analysis can be generalized when the spectrum satisfies 



(la.i) in 5 dimensions : ^ — ^ < 



31 



(F.2) 



which is the analogue in five dimensions of the first condition (3.13). One finds that the small 
terms of order (9(e^^") are again attracting the Universe towards phase 2. The difference 
with the case (I&.i) above is that along the RDS in phase 3, z is not stabilized but only frozen 
at any value such that e^^° ^ ij^ As before, we expect the system will enter phase 2 (as 
could be checked numerically), where ( freezes. However, the final space-time dimension in 
phase 2 can be four (with z stabilized) or five: 



V 



4 dimensions if < ^ — } 

15 + rirp 



< 



1 

'31 



5 dimensions if 



+ nip 



< 



1 

15 



(F.3) 



"'^^To observe the attraction from phase 3 to phase 2, we take the hmit e^^ ^ 1 in Eqs (5.17) and (5.19). 
One finds a dominant term 0{e~^~^) plus two subdominant monomials of order e^^'^^^ and e"*^, respectively. 
Neglecting these two terms, an RDS is found. However when these residual forces are taken into account, 



one finds = e^y"^^ ( 1 + =|±41 - e^^^f^ f 
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as a consequence of the dynamical decompactification of the direction 4 in the later case. 

• The above analysis can be applied similarly to case (II. i) when ny<0 {i.e. cases (Ila.i) 
and (II6.i)). The Universe re-enters phase 2 where ( freezes, z is stabilized, and the final 
RDS attraction is four-dimensional. 

• Let us apply the same techniques to case (I.ii). First, we consider models where 



(16. ii) in five dimensions : < — 

31 rin 



< 0, 



(F.4) 



The difference with the discussion of the case {Ib.i) is that we have a subdominant term 
in —T~'^J^ of order zero in e^, k{z) := n*j,frp[z) + nyfy{z). Moreover, the "residual force" 
induced by k, can either be positive or negative when Uy < 0, depending on the value of 
z (instead of being always positive). This could imply that in some models we may not 
be attracted back to phase 2. However, reasoning as in Sect. 3.3.3 of one finds that 
K induces a back-reaction on z and we expect that the residual force on y will always end 
by increasing it. In any case, this fact is always true in the explicit models considered in 
the present paper, since the condition (F.4) implies riy = and thus k(z) > for all z 



(see Eqs (2.12) and (2.15)). The Universe is attracted towards phase 2,-'^^ where it enters a 
four-dimensional RDS. 



Again, the analysis can be generalized when we have 

(la.ii) in five dimensions : 



< 



1 

31 



(F.5) 



As usual, y < increases and we approach phase 2. Once there, the final RDS is either 



four or five dimensional, as indicated in Eq. (F.3). Note that our analytic discussion does 
not cover the situations where riy + ny < 0, with riy > 0. However, we expect the attraction 
from phase 3 to 2 to take place as well. 



for K 



The discussion in case (H.ii) for riy < is as in case (I.ii), with a simpler expression 

fj. — 



n 



T 



1 
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The conclusions are identical, except that the final RDS in phase 2 



is always four dimensional. 



13 



Quantitatively, one finds e(^y) — 



when K is neglected in phase 3 
""^^Reasoning as in footnote 



_9_ raj/r+»V/v 



> 0, where is the value where z is stabilized 
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one finds e 



iv) 



c5zo 3 c 
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One finds in case (116. ii) e(j,) 



.v„ 9 »t(1-|)':4 
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5 (K 
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J_ 1 _i_ pjyo 183 £5 !hr ^ A 
31^^^ 80 ce ^ 



> 0, where z'^ is the value where z is sta- 
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